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UNITARY AND ORTHOGONAL EQUIVALENCE OF SETS OF 

MATRICES 

NAIHUAN JING 


Abstract. Two matrices A and B are called unitary (resp. orthogonal) equivalent 
if AU = VB for two unitary (resp. orthogonal) matrices U and V. Using trace iden¬ 
tities, criteria are given for simultaneous unitary, orthogonal or complex orthogonal 
equivalence between two sets of matrices. 

1. Introduction 

Let {Ai,A 2 , ■ ■ ■ , Ak} and {Bi, B 2 , • • • , B^} be two sets of complex matrices of size 
m X n. We say that {A} and {R} are (simultaneous) unitary equivalent if there exist 
two unitary matrices U, V of respective dimension such that 

(1.1) V*AiU = Bi, i = l,---,k, 

where * means the transpose and conjugation. Two sets of real (resp. complex) ma¬ 
trices {Ai}, {Bi} of the same size are said to be orthogonal (resp. complex orthogonal) 
equivalent if there are orthogonal (resp. complex orthogonal) matrices P and Q such 
that Q^AiP = Bi for all i. When the matrices are square matrices and U = V, then 
two sets are said to be (simultaneous) unitary similar. We also say that {Ai} are si¬ 
multaneous complex orthogonal similar to {R} if AiO = OBi for a complex matrix O 
such that 00^ = 0^0 = L 

Simultaneous unitary similarity has been an important problem in representation 
theory. When the set of matrices has no group structure, one of the hrst nontrivial 
results was given by Specht [12], and the question has been studied by many people, 
see in particular [1] for its history and difficulty. In [7| algorithms are given to related 
problems on simultaneous unitary similarity or congruence for complex square matrices 
(see also 0). Recently, important applications are found in quantum computation, 
where a key question of local unitary equivalence between two quantum states has been 
reduced to simultaneous orthogonal equivalence between two sets of real matrices (cf. 
|2|)- Current resurgent interest in quiver theory is also a reflection of the importance 
of the problem. 

Geometrically fll.ip represents the matrix representation of a set of linear transfor¬ 
mations under common change of orthonormal bases in the domain and range spaces. 
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If we relax fll.ip to only require that AiP = QBi for two non-singular matrices P, Q of 
respective dimensions, then Ai and Bi are two matrix representations of a set of linear 
transformations Lj under different bases. The special case of two linear transforma¬ 
tions {k = 2) was answered by Kronecker’s theory [6] of the matrix pencil Ai + XA 2 
using elementary divisors. In further special cases of symmetric and antisymmetric 
matrices there are canonical forms under similarity [13]. But for the corresponding 
unitary problem appropriate verihable conditions of P and Q are needed. In various 
applications in quantum computation, one is concerned with the problem of how to 
judge simultaneous orthogonal equivalence rather than hnding the intertwining matrix, 
as the hnal solution relies on the relevant problem in invariant theory. 

In this note we will give criteria for fll.ip and its analogues and show that two 
sets of complex (real/complex) matrices are unitary (orthogonal/complex orthogonal) 
equivalent if and only if the corresponding traces of words in AiA* (resp. are 

invariant. Usually the real version of a problem is harder, fortunately there is a satis¬ 
factory solution for almost all of our statements thanks to a simple argument to pass 
from the complex numbers to the real numbers. Our approach uses some results from 
semigroup theory. As an application, we also obtain an alternative version of Albert’s 
criterion for simultaneous similarity of symmetric matrices fRemark I2.6p . which avoids 
the complicated special cases in the original argument. 

We remark that a care is made to be self-contained and to indicate the generalization 
to semigroups from the corresponding results in hnite groups. 

2. Square matrices 

We hrst consider the unitary similarity (resp. orthogonal similarity) of two sets of 
complex (resp. real) square matrices. Let {Aj} be a set of complex (resp. real) matrices 
of same dimension. We say that {Ai] is hermitian closed (resp. transpose closed) if 
each A* (resp. A*) is contained in the span (Aj) of the set {Aj}. Two sets {Aj} and 
[Bi] of complex (resp. real) square matrices of equal dimension are said to be unitary 
similar (resp. orthogonal similar) if there exists a unitary (resp. orthogonal) matrix U 
such that U*AiU = Bi for all f, or equivalently U intertwines Ai and Bi for all i. 

Let iS be a set of matrices of equal dimension, a (product) word w{S) in the alphabet 
S' is a matrix product xy ■ ■ ■ z, where x,y,..., z are arbitrary matrices in S. We use 
W{S) = {tc(S)} to denote the set of words in the alphabet S. 

To give our hrst main result we need to recall some basic notions of semigroups and 
their modules. A set G is a semigroup if there is a closed binary operation on G that 
satishes associativity and has an identity. If every element of G is invertible, then 
G becomes a group. We follow [H I] to say that a vector space U is a G-module if 
there is an action of G on U: G x U 9 (x, u) 1 —)■ x.v G V such that x.{y.v) = {xy).v 
for all X, y G G,v G V. We assume that all modules considered in this paper are 
hnite dimensional left modules, but we do not assume that the semigroup G is hnite. 
Equivalently if G is a (semi)group of linear transformations on V, then U is a G-module 
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with the action given by the transformation. In this case, the entries of the matrix 
representation are called the coordinate functions of G. 

Two G-modnles Vi are called equivalent or isomorphic, denoted by Vi ~ V 2 , if there 
is a linear map (f : Vi — > V 2 snch that x.(j){v) = (f){x.v) for all x G G, n G V. In terms 
of the matrix representation, this means that xP = Px (Vx G G C End(y)) for some 
non-singnlar matrix P. 

The notions of reducible, completely reducible and decomposable G-modules can be 
dehned as in the sitnation of gronp modnles (cf. [HI I]). The most nsefnl ones for ns 
are (i) V is an irredncible G-modnle if V has no non-trivial G-snbmodules. It is known 
that the Schnr lemma holds in this case, i.e., the only G-homomorphism between two 
irredncible G-modnles are scalar homomorphisms. (ii) E is a completely redncible 
G-modnle ii V ~ Vi © • • • © 14- where 14 are irredncible. 

If G is a semigronp of linear transformations on V, then the associated G-module 
V is redncible iff nnder the matrix representation, every element of G is nniformly of 


the block triangnlar form x 


xi 0 
X 3 X2 


of hxed shape, where xi and X 3 are sqnare 


matrices of size smaller than that of x; and V is decomposable if all x 


Xi 

0 


0 

X2 


with 


the hxed shape. 

A G-modnle V is called unitary if there exists a positive-dehnite hermitian form ( , ) 
on V snch that (xn, xv) = {u, v) for any x G G, and arbitrary u,v & V. Eqnivalently the 
corresponding linear transformations of the elements of G are represented by nnitary 
matrices. In the case of a transpose-closed G-modnle, this is eqnivalent to the existence 
of a G-invariant symmetric bilinear form on V. 

The character x(E) of G-module V is the trace function x(E)(x) = trv{x), x G G. 
We will show that two irreducible G-modules are equivalent iff their characters are the 
same (Theorem 12.1[) . 

In hnite group theory, every complex module is completely reducible (Maschke’s 
theorem). For the semigroup theory, this is not true in general, but we still have some 
form of Maschke’s theorem for (possibly inhnite) semigroups. 

The following theorem was mainly due to Frobenius and Schnr [5] for linear trans¬ 
formations. We remark that the statements hold for semigroups of transpose closed 
linear transformations as well. 


Theorem 2.1. Let G be a semigroup. (1) If V is a unitary G-module, then V is 
completely reducible. (2) Two completely reducible G-modules are equivalent if and 
only if their characters are equal. 

Proof, (i) The idea of the proof is to show that if V is unitary reducible, then it is also 
unitary decomposable, i.e. the submodule has a complementary submodule or it is a 
direct summand. Using the matrix representation, this boils down to the fact that a 
block triangular matrix must be block diagonal if it is invariant under *-operation. 
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(ii) Since completely reducible modules are direct sums of irreducible modules, by 
adding necessary irreducible summands with possibly zero multiplicity, we can write 
the completely reducible modules as 

[/ = y®'>i ©...©y6, eZ+ 

where Vi are pairwise inequivalent irreducible G-modules and assume that xiY) = 
x{U). Taking characters we have 

(ai - bi)x{Vi) H-h (a^ - &r)x(K) = 0. 

By a theorem of Frobenius and Schur [5] (cf. [3l Th. 27.8]) the coordinate functions 
of pairwise inequivalent irreducible modules for semigroups are linearly independent, 
their characters are thus linearly independent. Therefore a* = for z = 1,... ,r, i.e. 
1 / ~ [/. □ 

We now come to our hrst result. It was announced in HU. The special case of a 
pair of matrices is the Specht criterion na for {A, A*} and Since the basic 

results and notions for semigroups are prepared above, a simple proof can be furnished 
as follows. 

Theorem 2.2. Let {Ai} and {Bi} be two sets of n x n hermitian closed complex ma¬ 
trices. Then {Ai} and {Bi} are unitary similar if and only if tr{w{{Ai}) = tr{w{{Bi}) 
for any word w. 

Proof. Let the index set of {Ai} be /, and consider the free semigroup G generated 
by 2 / 1 , •• • ,y\i\- The assignment yi i-A Ai (resp. Bi) defines a representations Va (resp. 
Vb) oi G on the n-dimensional space V = C^. Their associated (j-modules are also 
denoted by Va and Vb as well. 

The hermitian closedness implies that both modules Va and Vb are unitary com¬ 
pletely reducible modules for the free semigroup G. Since the trace of every word in G 
are the same, Va and Vb have the same character x(Va) = x(Vb). By Theorem 12.11 it 
follows that Vh — Vb, so there exists a non-singular matrix P such that 

( 2 . 1 ) AiP = PBu t = 

Since {A^} and {Bi} are hermitian closed, we also have 

(2.2) P*Ai = BiP\ z = l,---,fc. 

Then PP*Ai = PBiP* = AiPP*. As PP* is positive semi-definite, we let (PP*)V2 
be the square root of PP*. Then (PP*)i/2 also commutes with A* for all i. Let U 
be the unitary part of P in the polar decomposition. Subsequently U = (PP*)“V2p 
intertwines with A* and Bi for all z, i.e. AiU = UBi for i. □ 

It can be shown that only finitely many trace identities are needed. For a pair of 
matrices, the bound of the word length is nw + i + P!- For sets of real 
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symmetric matrices under Jordan-closeness, the bound can be improved (see Remark 
I2.6p . The real version of the theorem also holds. 

Corollary 2.3. Let {Ai} and {Bi} be two sets of real square matrices of the same 
size and assume that both sets are closed under transpose. Then {A} and {Bi} are 
orthogonal similar iff triw^Af)) = tr{w{Bi)) for any word w in respective alphabets, 
and iff there is a real matrix P such that PAiP~^ = Bi for all i. 

Proof. By Theorem 12.21 the trace identity implies there exists a unitary matrix U such 
that AiU = UBi for any i. Let P and Q be the real and imaginary part of U, then 
one of them is nonsingular, say P. Taking the real part of the intertwining equation, 
we obtain that AiP = PBi for any i. Let O be the orthogonal part of PP^ in its polar 
decomposition, then the same argument of Theorem 12.21 gives that AiO = OBi for all 
i. □ 


If we replace the hermitian closeness condition by transpose closeness, then we have 
the following result. 

Theorem 2.4. Let {Af} and {Bf} be two sets ofnxn transpose closed complex matri¬ 
ces. Then {Af} and {Bi} are complex orthogonal similar if and only if tr{w{{Ai}) = 
tr{w{{Bi}) for any word w in respective alphabets. 

Proof. We noted that Theorem 12.21 holds for semigroups of matrices that are transpose 
closed. The same proof of Corollary 12.31 can be repeated to derive the result. □ 


Remark 2.5. The criteria do not hold without the condition of hermitian or transpose 


closeness. For example let A = 
not similar to B. 


1 1 
0 1 


,B = 


1 0 
0 1 


. Then tr{A'^) = tr{B^), but A is 


The following application modihes Albert’s criterion on simultaneous similarity of 
real symmetric matrices. 

Remark 2.6. Let {A} o,nd {A} be two sets ofk real symmetric matrices. Suppose that 
the traces of any corresponding words in A and Bi are identical. By a classical fact that 
any symmetric polynomial is a polynomial in the power-sum symmetric polynomials, 
it follows that det{xl — A) = det(x/ — B) for any linear combination A = ^iAi 
and B = Yl!i=i^iAi, where x,Xi are indeterminates. Albert’s criterion [ 1 ] says that 
{Ai} is “almost” simultaneous orthogonal similar to {Bf} if the determinant identity 
holds plus that {A} and {Bf} are Jordan-closed, i.e. closed under anti-commutators 
{X,Y} = XY -\-YX for any members of the sets. Here “almost” means that there 
are some counterexamples of degree 2 simple Jordan algebras of dimension Aqi + 2. 
Moreover, [1] showed that these “counterexamples” can be removed if certain monomials 
in {A} and {A} of length Aqi + 2 have the same traces. In particular, one only needs 
to verify the trace identity for word length < max{n,4:qi-\-2) to ensure the simultaneous 
similarity of {A} to {A} under the condition of Jordan-closeness. 
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3. Rectangular matrices 

The following elementary fact is needed for fnrther discussion. 

Lemma 3.1. Let a,b (resp. be two column vectors. Then aa* = bb* iff there 
is a complex number 9, \6\ = 1 (resp. 6 = ±1) such that a = 9b. 

Proof. Let a = (oi, • • • , OnY and b = ( 6 i, • • • , Clearly a* = 0 iff fej 
can assume that Oi, • • • , are non-zero numbers, and a^+i = ■ ■ • = 
hip = hip h 0) ^6 can write a* = 9ibi with \9i\ = 1 for each i = 1, 
aiOj = bibj h 0 imply that 9i9j = 1 for any i,j = 1, - ■ ■ k. Therefore 9i = ■ 

Then we have that a = 9b. 

We have the following generalization. 

Lemma 3.2. Let A, B be two m x n-matrices over C (resp. W). Then AA* = BB* 
iff there is a unitary (resp. orthogonal) matrix V such that A = BV. 

Proof. The sufficient direction is clear. On the other hand, let U = hi,-- - ,Um\ be 
the unitary matrix of a basis of orthonormal eigenvectors of AA* = BB*. Suppose 
the hrst r eigenvectors have nonzero eigenvalues af, then Vi = a~^A*Ui exhaust all 
eigenvectors of A*A with non-zero eigenvalues by the singular value decomposition 
(same eigenvalue af). Extend {ui, - - - , Ur} into a unitary matrix Vi = hi)''' > ^n] of 
a basis of orthonormal eigenvectors for A*A. Then AVi = UD,A*U = ViD^, where 

D = Q j and Di = diag{ai, - - - jCr^). Similarly there exists another unitary 

matrix V 2 such that BV 2 = UD,B*U = V 2 DK Subsequently A = D 

Theorem 3.3. Let {Aff and {Bf) he two sets of complex (resp. real) m x n-matrices. 
The following are eguivalent. 

(a) the set {At} is unitary eguivalent (resp. orthogonal eguivalent) to the set {R*}; 

(b) the set {AiA*\i < j} is unitary similar (resp. orthogonal similar) to the set 
{B,B*\i<j}; 

(c) tr w{{AiA*}) = tr w{{BiB*}) for any word w in respective alphabets. 

Proof. Equivalence of (b) and (c). Suppose (b) holds, then there is a unitary matrix 
Lf such that U*AiA*U = BiB* for any i < j. Taking * we also have Lf*AjA*Lf = BjB* 
for any i < j, then Lf*AiA*U = B^B* hold for any i,j. Therefore the set {AiA*\i,j} 
is simultaneous unitary equivalent to the set {RjR*hi}, so (c) holds by Theorem 12.21 
The converse direction is guaranteed by Theorem 12.21 

(b) clearly follows from (a). We now show that (b) implies (a). Let A be the block 
matrix dehned by A* = [Al, -- - ,Al]. Then the block matrix AA* = [A^j] with 
(f,j)-ccitry being an m x m matrix AiA*. By the argument above (b) implies that 
there is a unitary matrix U such that 

(3.3) UAiA*U* = BiB* 


= 0 , so we 
= 0. Since 
- - - ,k. But 
■ ■ = 9k = 9. 
□ 
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for any i,j. Then 



-UA{ 



-B{ 

(3.4) 


[A\U* ■ 

• KU*] = 



UAk. 



-Bk. 


matrix V such that 


-UA{ 


1- 

cq • 

1_ 

UAk. 


.Bk. 


V. 


[b; Bi] 



ruAi] 


[Ba 

Using Lemma |3.2| for the block matrices 


and 



[uAk\ 


.Bk. 


, we get a unitary m x m- 


Subsequently UAi = BiV for all i. 

The real case follows by a similar argument in view of Corollary 12.31 □ 


The following result is clear from our discussion. 


Theorem 3.4. Let {A} cmd {B^} be two sets of complex mxn-matrices. The following 
are equivalent. 

(a) the set {A} is complex orthogonal equivalent to {Bi}; 

(b) the set {AiAj\i < j} is complex orthogonal similar to {BiBj\i < j}; 

(c) tr w{{AiA^-}) = tr w{{BiB}}) for any word w in respective alphabets. 


ACKNOWLEDGMENTS 

The author thanks Professors M. Putcha and V. V. Sergeichuk for interesting discus¬ 
sions. He is also indebted to the referee for suggestions that have improved the paper. 


References 

[1] A. A. Albert, On the orthogonal equivalence of sets of real symmetric matrices, J. Math. Mech. 
7 (1958), 219-235. 

[2] Y. A. Alpin, K. D. Ikramov, On the unitary similarity of matrix families, Mat. Zametki 74 
(2003), no. 6, 815-826; transl. in Math. Notes 74 (2003), no. 5-6, 772-782. 

[3] C. W. Curtis, 1. Reiner, Representation theory of finite groups and associative algebras, John 
Wiley & Sons, Inc. New York, 1962. 

[4] S. Friedland, Simultaneous similarity of matrices, Adv. Math. 50 (1983), 189-265. 

[5] G. Frobenius, I. Schur, Uber die Aquivalenz der Gruppen linearer Substitutionen, Sitzungsber. 
Akad. Berlin (1906), 209-217. 

[6] F. R. Gantmacher, The theory of matrices I, II, Translated from the Russian by K. A. Hirsch. 
Reprint of the 1959 translation. AMS Chelsea Publishing, Providence, RI, 1998. 

[7] T. G. Gerasimova, R. A. Horn, V. V. Sergeichuk, Simultaneous unitary equivalence. Linear 
Algebra Appl. 438 (2013), 3829-3835. 

[8] N. Jacobson, Basic Algebra I, Second ed. W. H. Freeman and Company, New York, 1985; 
Basic Algebra II, 1989. 


















NAIHUAN JING 


[9] N. Jing, M. Li, X. Li-Jost, T. Zhang, S. Fei, SLOCC invariants for multi-partite mixed states, 
J. Phys. A 47 (2014), No. 21, 215303. 

[10] C. J. Pappacena, An upper bound for the length of a finite-dimensional algebra, J. Algebra 197 
(1997), 535-545. 

[11] N. A. Wiegmann, Necessary and sufficien conditions for unitary similarity, J. Austral. Math. 
Soc. 2 (1961/1962), 122-126. 

[12] W. Specht, Zur Theorie der Matrizen II, Jber. Deutschen Math. Verein. 50 (1940), 19-23. 

[13] R. C. Thompson, Pencils of complex and real symmetric and skew matrices. Linear Algebra 
Appl. 147 (1991), 323-371. 

Department of Mathematics, North Carolina State University, Raleigh, NC 27695- 
8205, USA 

E-mail address: jing@math.ncsu.edu 


